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1. Introduction 



Chronoamperometry, where a current versus time 
curve is measured after a step potential is applied to the 
electrode, is one of the important techniques in electro- 
analytical chemistry. The analysis of chronoamperome- 
try at a shrouded circular disk electrode (see Fig. 1) 
involves solving boundary value problems with 
boundary and initial conditions that give rise to discon- 
tinuities in concentration and flux distribution at the 
dislc/shroud boundary. Discontinuous integrals of Bessel 
functions, which are solutions to the diffusion equa- 
tions, have this property; some of them have been used 
in heat conduction in cylindrical geometry [1]. They 
have also been nicely applied in treating steady state 
problems in electroanalytical chemistry [2,3]. However 
applications of these integrals in transient processes, 
such as chronoamperometry, at the shrouded disk elec- 
trode become more difficult and may involve current 
functions in the Laplace domain that can not be easily 



inverted analytically [4,5]. We shall report in this paper 
some of the mathematical results we obtained in apply- 
ing a class of these discontinuous Bessel integrals to a 
boundary value problem related to that of chronoamper- 
ometry and compare our results with some of the results 
in chronoamperometry. 



2. Concentration Distribution 

Consider a disk electrode with radius a, and let the 
origin of a cylindrical coordinate system be placed at 
the center of the disk with z-axis perpendicular to the 
disk. Also consider a reaction at the electrode, 
Ox+ne^Red, and let the concentration of Ox, 
C{r,z,t), be C° throughout the solution initially. At 
time t=0, let the potential be stepped to a large negative 
value such that the concentration of Ox at the electrode 
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Fig. 1. Shrouded electrode. 

is zero. Then C(r,z,t) satisfies the diffusion equation, 
dC(r,z,t)/dt=D[dVdr^+r~'d/dr+dydz^]C(r,z,t) 

with the following initial and boundary conditions: 
C(r,z,0)=C° att=0 
C(r ,0, 0=0 for 0<r<a , f>0 
\\mC{r,z,t)^C° 

limC(r,z,t)^C° 

and D denotes the diffusion coefficient. For the com- 
monly used shrouded electrode there is an additional 
boundary condition that the flux into the shroud is zero. 
The mixed boundary conditions of zero flux at the 
shroud and zero concentration at the disk make a more 
difficult mathematical problem. Various techniques 
have been applied to this boundary value problem with 
various degree of rigor [6-9] . We shall consider the lim- 
iting case where the thickness of the shroud approaches 
zero and ignore the zero flux condition for the moment. 
Let us define the following dimensionless variables: 
x=r/a,y=z/a, u=Dtla^, and c{x,y,u)= C{r,z,t)IC°—\; 
then the above diffusion equation and the boundary con- 
ditions become, 

dc(x,y,u )/du=(dydx^+x~ ' d/dx+ d^/dy ^)c(x,y,u) ( 1 ) 

c(x,y,0)=0 (2) 



c(x,0,m)=-1 forO<x<l 
lime (x,y,u )=limc (x,y,u)=0 



(3) 
(4) 



Taking Laplace transforms of Eq. (1) with the initial 
condition Eq. (2), we obtain, 



(aVax'+x" ' d/dx+d^/dy'--p )cix,y,p)=0 



(5) 



where p denotes the Laplace transform variable. Equa- 
tion (5) is solved_by the method of separation of vari- 
ables by letting c(x,y,p)=X(x,p)Y(y,p). Substituting 
this into Eq. (5) the following two ordinary differential 
equations are obtained 



(d'/dx'+x-'d/dx+f)Z(x,/7)=0 
(dVdy^-X^)Y(y,p)=0 



(6) 
(7) 



where ^^ is the separation parameter and \^=^^+p. The 
solutions to Eqs. (6) and (7), with the boundary condi- 
tions of Eq. (4), are the Bessel functions Joi^x) and the 
exponentials exp(— Aj), foT X(x,p) and Y(y,p) respec- 
tively; we denote the set of product solutions to Eq. (5) 
by g (^) Jo(^x)e.xp(— \y), and try to satisfy the discon- 
tinuous boundary condition Eq. (3) by combinations of 
product solutions with various ^. To satisfy Eq. (3), we 
try to make use of some of the discontinuous integrals 
of Bessel functions [10] such that J g (^) Jo(^x)d^— \lp 
for 0<x<l. One of the well known discontinuous Bessel 
integral with ^(^)=A:sin^/^ will give a constant for 
0<x<l, but it does not give a value of for all x>l at 
t=0, Eq. (2). We shall try a class of discontinuous Bessel 



integrals of 



Weber with 
as an 



g(^)=^^(a/P"', writ- 



ing c{x,y,p) as an integral, i.e., c{x,y,p)= 
^J^'"%(^)-/o(^x)exp(- Aj)d^. At j=0. 



c(. 



xAp)=k^ 



^'-%(a/o(^x)d^ 



=k{{\ -x^yiy-'IFiix) for 0<x<l 



=0 



for l<x<oo 



(8) 



For iJi=\ and k=—\lp, the boundary condition at the 
disk, Eq. (3), is satisfied and also c=0 for x>l. This is 
the only member of this group of functions that satisfies 
the condition that c(x,0, is uniformly equal to —1 at 
all X less than 1. However if instead we only require that 

Eq. (3) be satisfied on the average, i.e., 

I I 

(c(x,0,M))=/c(x,0,M)xdx//xdx=— 1, then the other 
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members of this group satisfy this condition with 
k=—2'^~^r(ix+l)/p. We shall proceed with this and take 
the inverse transform to obtain, 

c(x,y,u)=L-\-2>'-T(fji+l)/p) 
je~''UOJoax)exp[-y(p+ey''d^ 



,+1) r ^(g)/o(gx 



Hm+D I 



+e-'^erfc(y/2V^-V^^)]d^ (9) 

3. Flux Distribution 

The flux at y=0 is given by 



J t J 



dyr°~ 2 
1 



'TTM 



For ease of calculations, the first integral is written as 
? ^(g)/o(gx>r/(\4g) , f ^(g)/o(gx) , 

I pZ^^2 tlb= I >;fl"2 tlfc 



p^(g)/o(gx>r/c(\4g). 



t;U,-2 



-d^ 







2"-T(3/2) 



■/^ 



■F(3/2,3/2-^;l;x") 
d^ x<l 



r(3/2-^) 

'/^(g)/o(gx>r/c(\4g) 



tii-i 



2-ix 



■F(3/2,3/2;^+l;x-") 



-f 



xT(^+l) 



(11) 



for ^(,>1, where F(a,i';c;z) denotes the Gauss hyperge- 
ometric series [11]. Figure 2 illustrates (dc/dy) at ^=0 
as functions of x and m, for ix=2 and 3. For x>l, (dc/ 
dy)y=o is nearly zero, this is particularly true for small u 



and forx away from the neighborhood of x=l. And as u 
approaches 0, the first integral in Eq. (10) approaches 

and the second integral tends to /^'~''A(^)/o(^x)d^ 



[12] which is identically forx>l, Eq. (8); and the sec- 
ond term in Eq. (10) numerically becomes almost for 
x>l as M^O. Thus the concentration distribution given 
by Eq. (9) approximately satisfies the condition of zero 
flux at x>l, i.e., at the shroud. There is a singularity in 
(dc/dy)y=o at the point x=l for the case of ix=2, because 
F(a,b,c;l) is divergent for /U-^2; nevertheless, the cur- 
rent, which is the integral of the flux over the disk area, 
is regular for ix=2. 

4. Current as Function of Time 

The current is given by i(t)=nFD j(dC /dz)z=o'2'nrdr , 
integrating over the disk, where F denotes the Faraday 
constant. Defining a dimensionless current function 
<p(u)=i(t)/2'TmFDC°a we get 



1 



4>(u)= (dc/dy)y=oxdx 



2' '>" 



^(g)/o(gx)er/(g\4) , 



1 



1 



? ^(g)/o(gx)e-"^ 



} 



■d^ 



xdx 



^(g)/i(g)gr/(g\4) 



d^ 



V TVUJ 



(g)/i(g)e-"^ 



e 



■d^ 



(12) 



When M— >0, i.e., when t«a^lD, the second term on the 
right hand side of Eq. (12) dominates and 

limV^<^(«)=i?^f4(|lM)d^ 1 

"^0 ^' 2'"''V^J ^'^ 2V^ 



or, in terms of the physical variables, 
nFC°Vid5a^ 



Dt 



lim i{t)= 
^0 



Vt 



which is the Cottrell equation. 
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Fig. 2. Fluxes as functions of radial distance and time. From top to bottom: Vm=0.05, 0.1, 0.5, 1.0. 

^(V^)Ji(V^) 



Series expansion for large u: The last integral in Eq 
(12) can be evaluated for large u from the following re- 
lation [13], 



<Pi{ 



»4i 



(g)./i(g)e 



-ue 



e 



■d^ 



-2-S 



-(m+l) 



rifi+2m+2)i-4u) 
;^ r(fjb+m+2)r(fjb+m+l)r(m+2) 



and the first integral on the right hand side of Eq. (12) 
from 

^.(«)= f-^^^^-^g!f^^^ dg 



2-T(fji-l) 

~r(^+i/2)r(^-i/2) 

Then the current, Eq. (12), is evaluated by 



u 

-F= —F=id(p2ia)/da)da. 



4>{u)= 






niJi+1) 



Am-1) 



r(ix+-)r(ix--) 



riiJi+2m)i-4u)' 



(2m+l)r(fjb+m+l)r(fjb+m)m ! 



which is good for large value of u, i.e., for long time. 

Asymptotic expansion at short time: We now seek an 
asymptotic expansion at small ufoT (j)(u). Let us change 
the integration variable in (p2{u), the last integral in Eq. 
(12), from^t0 7=^^ 



<p.(M)=l/2j^^ 



,(H+l)/2 



e "Mt. 



Thus (p2{u) may be regarded as the Laplace transform 
of the function 

The MelUn transform of/(T) is [14], 



(i/2)-^^-^''-^^r(^-2z+2)r(z) 



r(^-z+i)r(^-z+2)r(2-z) 

which is absolutely convergent for 0<Re z< (/a+3)/2. We 
therefore apply the results of Handelsman and Lew for 
the asymptotic expansion of Laplace transform near the 
origin [15,16] to obtain 

c+ico 

^2(M)=(i/)(M)=(2T70"' J F{z)r{\-z)u'-'dz 
0<c<l 

(-2)""^^ 



2"" 



hE 



^ li+m+2 . (m+„-i)/2 



r(u+l) fl ,„,u-OT+2, „,u-m, , , 



and the first integral, ^i(m), is derived from 

u 

(Pi(u)= ^ -^(d^2(a)/da)dQ;. 
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Therefore, at small u the current, Eq. (12), is given by. 



<A(m)= 



2V^ 



-2''r(^+i) 



1 



i-iTrc 



ix,+m+l (m+„-i)/2 



m !r(^^^^)r(^)(^+m+l)(^+m) 



Several <^(m) versus 1/ vm plots are shown in Fig. 3, 
including the Cottrell equation for the one dimensional 
diffusion, which is the straight line through the origin. 
The intercepts of the curves at the vertical axis show the 
steady state current, (j)(u—^°°), for iJi=2 and 3, and they 
have the values 0.848 and 0.678, respectively; while 
from numerical simulation of the chronoamperometry 
at inlaid disk electrode [17,9], <^(m->co)=0.637 which 



corresponds to a steady state current of 
i(t-^oo)=4nFDC°a. As fjb increases (j){u-^°°)= 
r(iJi+l)r(iJi—l)/2r(iJi+j)r(iJi—j) decreases monotoni- 
cally to 0, i.e., cj)(u) approaches Cottrell behavior as /x 
becomes large. In Fig. 4, yu(j)(u) for (j^2 and 3 are 
plotted versus vm to magnify results at small u. The 
results from numerical simulations for shrouded disk 
electrode are also included, curve b. Discontinue Bessel 
integral was used differently in chronoamperometry by 
Fleischman and Pons [5], their result is also shown in 
Fig. 4, curve c. 
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Fig. 3. (^(m) as functions of 1/V«. 



Fig. 4. \/u<j>(u) as functions of V«. (a) yLi=2, (b) digital simulation, 
(c) Ref. [2], (d) /Li=3. Horizontal line is Cottrell one dimensional 
diffu- 
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